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We study the problem of witnessing entanglement among indistinguishable particles. For this 
purpose, we derive a set of equations which results in necessary and sufficient conditions for probing 
multipartite entanglement between arbitrary systems of Bosons or Fermions. The solution of these 
equations yields the construction of optimal entanglement witnesses for partial and full entanglement 
in discrete and continuous variable systems. Our approach unifies the verification of entanglement for 
distinguishable and indistinguishable particles. We provide general solutions for certain observables 
to study quantum entanglement in systems with different quantum statistics in noisy environments. 
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I. INTRODUCTION 

Nonlocal correlations among many particles or quan¬ 
tized fields are one key element of the quantum nature of 
physics m- Applications in metrology use this quan¬ 
tum feature to beat classical limitations [SHE]. Quan¬ 
tum entanglement has been also studied as a resource for 
quantum information technologies [9]. For a fundamen¬ 
tal characterization, a lot of attention has been devoted 
to verify entanglement between distinguishable particles 
(DP) or multiple degrees of freedom [lDl ITT] . 

Indistinguishable particles (IP), on the other hand, 
are indispensable for understanding the properties of 
many-particle quantum systems. For IP systems having 
a certain spin statistics |12j . however, even the notion 
of entanglement itself has no generally accepted defini¬ 
tion HSj. For example, the two-Fermion or two-Boson 
state, 

lt)®U>-U>® I t) =| t) A | 4-), m 

11) ® U> + U) ® 11) =11) VI ;>, u 

is a Bell-like entangled state using the tensor product (g), 
and, at the same time, it is a product state in the notion 
of the antisymmetric product A or symmetric product V, 
respectively. This ambiguity originates from the fact that 
the (anti)symmetrization requirement of (Fermion)Boson 
systems has formally the same structure as a nonlocal 
superposition. The left-hand side of Eq. |T]) is closely re¬ 
lated to the well established theory of entanglement be¬ 
tween distinguishable subsystems nn. Hence, we will fo¬ 
cus on the right-hand side pSHE]. That is, for Fermions 
and Bosons, the states 0 are separable product states in 
the exterior algebra and symmetric algebra, respectively. 
In general, we will focus on the following question: ”How 
does one certify quantum entanglement which does not 
rely on (anti)symmetrization?” 

For bipartite pure states, the relation between entan¬ 
glement for DP and the tensor product is represented by 
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the Schmidt decomposition 0. Whenever a single ten¬ 
sor product state is sufficient to expand a pure state, it is 
separable. In analogy, entanglement between IP is char¬ 
acterized by the Slater decomposition HSHIZ]. A pure 
Fermion or Boson state is separable, if it is a single an¬ 
tisymmetric or symmetric product state, given in terms 
of Slater determinants or permanents, respectively; see 
also HMU for a slightly different approach. A classical 
mixture of product states extends the corresponding defi¬ 
nitions to mixed quantum states. If such a representation 
is impossible, the state under study is entangled. 

Irregardless of the product for constructing compound 
Hilbert spaces, ®,A,V, one can detect quantum correla¬ 
tion via so-called entanglement witnesses [SEolEl], a 
witness is an observable which is non-negative for all sep¬ 
arable states, and may be negative for entangled states. 
Such criteria have been successfully applied to exper¬ 
imentally probe quantum correlations [22H27] , In the 
same dimensions, the characterization and application 
of entanglement in systems of IP gained an increasing 
importance during the past years; see, e.g., (28P401 . In 
particular, the entanglement of multiple qubits, realized 
in spin systems, has been investigated [HHS], and the 
question of how to extract entanglement for applications 
from systems subjected to the Pauli principle has been 
addressed [¥5J 0B] . Recently, a method for the construc¬ 
tion of optimized multipartite entanglement witnesses for 
DP has been proposed m and applied to perform a full 
entanglement analysis of experimentally generated mul¬ 
timode states IIE1- 

In the present contribution, we derive equations which 
allow the construction of optimized, necessary, and suf¬ 
ficient entanglement probes for Boson and Fermion sys¬ 
tems. The formalism is applicable to arbitrary numbers 
of particles, partially and fully entangled states, and dis¬ 
crete and continuous variable quantum systems. Our 
method unifies the detection of entanglement for DP 
and IP. Furthermore, we explicitly construct witnesses 
to demonstrate the strength of this technique to certify 
entanglement for different spin statistics in noisy envi¬ 
ronments. 

This work is structured as follows. In Sec. [TT] we dis¬ 
cuss the notion of multipartite separable states of IP. The 
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construction of corresponding entanglement witnesses is 
derived in Sec. |III| Bipartite examples are studied for IP 
and related to the systems of DP in Sec. m Section |V| is 
devoted to establish spin statistics independent witnesses 


in multimode systems. We conclude in Sec. VI 


II. A-SEPARABLE FERMION AND BOSON 
STATES 


The formulation of multipartite entanglement for DP 
is based on the tensor product structure of compound 
Hilbert spaces r H ® N . A TV-partite quantum state a is 
fully separable, if it can be written as a convex combina¬ 
tion of product states of the subsystems i 14, 


" f A Dt \ ’ * * ‘ P 

a= / dP(ai,...,a N )- - 

J (ai,..., 


\ai, .. ., a N ){ai,...,a N \ 
ciAr|ai, • • •, oat) 


( 2 ) 


Here, |ai,..., ajv) = |ai) ® ■ ■ • ® |ajv) are, in general, 
unnormalized TV-partite product vectors, and P is a clas¬ 
sical probability distribution. 

A fundamental postulate of quantum mechanics for 
Bosons or Fermions is that the quantum states are sym¬ 
metric or antisymmetric upon exchange of the subsys¬ 
tems, respectively. This restricts the physical states to 
the (anti)symmetric subspace of the TV-fold tensor prod¬ 
uct Hilbert space, 7T AAr , T~C /N C P® N . A projection from 
the tensor product space to these subspaces is given by 
the permutation operators II , 


n d 


I \ _ W (il)' 0 '' 

l a i, ■ ■ •, aN ) - 

0~£S N 


t(1) ) - * • , @cr(N))> (^) 


where \a\ and (±1)I CT I denote the parity and the sign of 
the permutation a £ Sn , respectively. Other projections 
might be similarly studied, which allow a generalization 
to other parastatistics. Now, the (anti)symmetric prod¬ 
uct states can be identified as [15] 


|ai) A ■ A |ajv) = n |ai,...,ajv), (4) 

|ai) V • • • V |ajv) = n + |ai,... ,ajv). (5) 

In general, for every vector | i/j) £ 'H® N , we get the 
(anti)symmetric vector in the projected subspace as 

IV^} = n ± |v>). (6) 

See Appendix A] for the symmetrization of operators. 

Equations ( 4| and © define fully or TV-separable 
Fermions and Bosons, respectively. More involved is the 
notion of AT-separable states; see mm for introduc¬ 
tions. In systems of DP a iV-separable vector \i/jx) £ Aa 
is defined as a product vector 

| 4’k) = |&i) ® • • ■ <S> | b K ) = |£»i, ■ ■ •, bx), (7) 

for positive integers (TVi,... ,Nk), \bk) £ H® Nk , and 
Y^k=i = N. The tuple (TV 1; ..., N K ) defines a par¬ 
titioning of the TV-fold Hilbert space. Further note that 


| bk) is, in general, not a product state in 'H® Nk . Finally, 
a A^-separable (anti)symmetric vector is defined as 

W’k) =n ± |^) =fl ± \b 1 ,...,b K ). (8) 

Since the permutation is applied in ([8]), the initial or¬ 
dering of the Hilbert spaces in (JT]) does not play a role, 

1. e., the partition (TVi,..., TV#-) relates to the Hilbert 
space sequence TT*^ 1 g> • • • g> 'H® Nk for DP or replacing 
(g) by A, V - for IP. Let us point out that (TV 1; ..., N K ) 
and (TV{,... ,N' K ) define, in general, different partitions 
if these tuples are not identical up to a permutation of 
indices. For example, the three partition (2, 3,1) of a six 
mode system describes the same partitioning as (1,2,3), 
but it differs from (2,2,2). 

An explicit example of partial separability of DP is the 
tripartite state: 

I#) = |0) ® (|1> g> |2) + |3) ® |4)) £ A 03 , (9) 

using a single-mode orthonormal basis {|0),..., |4)}. 
Similarly, one can construct partially separable states of 
DP in (§, e.g., 

|4/+)-|0)V(|1)V|2) + |3)V|4))£A v3 
and Id'”) =|0) A (|1) A |2) + |3) A |4)) £ FT A3 . 

In Appendix [B] we prove - independently from the 
method to be developed later on - that these projected 
states, |4' ± )=n ± |'P), are indeed partially separable, K = 

2, and not fully separable, K ^ 3. 

Based on the (anti)symmetric product, one gets a gen¬ 
eral definition of TV-separability for IP [33l [49]. Addi¬ 
tionally, a TV-Fermion or TV-Boson quantum state a is I\ 
separable, 1 < K < TV, if it can be written as a convex 
combination of (anti)symmetric product states, 



(bi ,..., bx |n=*= |6i, ..., bx) 

( 11 ) 


Exchanging the tensor product ® by either the symmetric 
product V or the antisymmetric product A, cf. Eqs. Q 
and <§> keeps the separability definition © of DP struc¬ 
turally preserved. If a state cannot be written according 
to definition (JTTJ) for K = TV, then entanglement be¬ 
tween IP is certified beyond any correlation that can arise 
from the (anti)symmetrization requirement of the quan¬ 
tum statistics itself. 


III. CONSTRUCTION OF ENTANGLEMENT 
WITNESSES 

A. Separation of entangled states 


Since Eq. (Ill defines a closed, convex set of states, 


the Hahn-Banach separation theorem is applicable 133 
EH- It states that for any closed, convex subset C of a 
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Banach space and a point, a; ^ C, there exists a linear 
and continuous functional / that separates these sets: 
f(x) > sup ceC /(c). In our case the Banach space is 
the set of Hermitian trace-class operators and the dual 
space is isomorphic to the set of bounded operators, i.e., 
f(p ) = tr (pL) for some Hermitian L. This means for the 
problem under study that for any iv-entangled state of 
IP, g = ]Q =t , exists a bounded Hermitian operator L, 
such that 

tr (Lg) > sup{tr(Lij) : for all a in Eq. ©}• (12) 

Related approaches to identify entanglement of DP can 
be additionally found in Refs. [20] HU 07] [55]. Since 
the separation theorem ensures the existence of such an 
operator, we have a necessary and sufficient condition 
in terms of observables L probing multipartite entangle¬ 
ment between IP in finite and infinite dimensional spaces. 
However, finding the proper L for a given state g and de¬ 
termining the least upper bound of the right-hand side 
of inequality © are cumbersome problems. In the fol¬ 
lowing, we will propose a method to address the latter 
aspect. 


B. Witness construction through optimization 


We define the abbreviation X for a Hermitian operator 
X. acting on 'hL® N , by the relation 

(x\X^\y) = (bi,, bk~i,x, 6/c+i, • ■ •, b K \X (17) 

x |&i,-•• A-i, 2 /) b k+ i,...,b K ), 


for all \x),\y) € V® Nk . Now, the derivative reads as 
dg d (&fc|(n ± Tn ± ) 5 -| 6 fe) 


0 = 


d(b k \ d(b k | (& fe |(n ± y 6 fe ) 

(tPLiPy^h) <6 fc | (n ± z,ft ± )^| fefc) 


(18) 


(n ± )d^>- 


(hW ± )^|ft fc ) < 6 *l(H ± )^| 6 fc ) 2 

Further, inserting the definition of g in Eq. ( |15[ ) yields 


0 = 




(19) 


Let us summarize the conducted optimization which 
results in Eq. (19). The derived set of algebraic equa¬ 
tions, 


(n ± M ± ) lr |6 i ) = &,•} for j = 1,..., K, (20) 


It is sufficient to take the least upper bound on the 
right-hand side of inequality ( 12 ) over all product vec¬ 


tors, being the extremal points of the given convex set 
of AT-separable states. Moreover, the operator II ± plays 
the role of the identity in the (anti)symmetric subspace. 
Combining these facts allows us to write condition (12) 


in terms of the expectation value of the entanglement 
witness operator: 

w± = Gn* -n ± M ± , (13) 

(bi, ■ ■ •, 6j{-|n ± Ln ± |6i, ...,b K ) 


G = sup • 


(&!,..., bx In* [&!,..., bx) 


, ( 14 ) 


where the least upper bound is taken over all product 
vectors ([T]); see also JT5, [TS]. Interestingly, the simple 
modification ^ —> 1 in Eqs. (13) and (141 yields the 


corresponding construction of witnesses for DP 07] [53]. 


The least upper bound in Eq. (14) defines an optimiza¬ 
tion of the Rayleigh quotient, 


_ (bi, ...,&AHn ± i(n ± |&i, ...,b K ) G 
9 ~ (b 1 ,...,b K \ii±\b 1 ,...,b K ) 


(15) 


define the separability eigenvalue (SEvalue) equations for 
IP. The common eigenvalue g is denoted as the SEvalue, 
and the (anti)symmetric product vector H ± \bi,... ,bx) 
is the corresponding separability eigenvector (SEvector) 
for IP. 

The SEvalue equations for IP represent a system of K 
coupled eigenvalue equations. Remarkably, it turns out 
that they have the same structure as the corresponding 
equations for DP m ■ For DP, we use the iV-fold identity 


1 that replaces the projector II ± in (20). Even though 


we find a strong relation to the case of DP, cf. mm , we 
want to point out that neither the distinguishable case 
includes the indistinguishable one, nor vice versa. This 
is due to the properties of the noninvertible operator n 111 . 

We found that the SEvalue g corresponds to an optimal 
expectation value of L for I\ -separable Boson or Fermion 
states. Therefore, we get the bound in the entanglement 
criterion ( 12 ) as 


sup{tr(L<x) : for all a in Eq. ©}=sup{.g}, ( 21 ) 

i.e., the initial convex optimization problem is solved by 
the largest SEvalue g of the equations (20). Now, the 


entanglement condition for the state g may be written as 


under the constraint that the denominator exists, i.e., 
II =t: |&i,..., bx) 7 ^ 0. The optimization is carried out as a 
derivative of the Rayleigh quotient: 


0 = 


dg 

Wi\ 


for 3 = 1,.... I\. 


(16) 


(L) = tr (gL) > sup{g}. (22) 


Alternatively, this condition can be written in terms of 
witnesses constructed from Eqs. |l3|) and (fl4|: 


W± = GL^ — n ± Ln ± , with G = sup{g} , (23) 
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which reads as tr {gW±) < 0. Similarly, a witness can 
be constructed using the lower bound of the Rayleigh 
quotient (1T5| as 


W± = IF 


L — inf {< 7 } 1 ft 


(24) 


IV. BIPARTITE EXAMPLE 

In a first application of our introduced method, we aim 
at witnessing bipartite entanglement. As our observable 
we may choose the rank one operator: 


Hence, by solving the algebraic problem (20) of ob¬ 
servables L, we are able to construct, in principle, any 
optimal entanglement witnesses for multiple correlated 
Bosons or Fermions. Moreover, the SEvalue equations 
for IP might be also used for a numerical optimization if 
an analytical solution is not available. Since the criterion 
in ( 22 ) and the witnessing approach are equivalent, we 
study from now on solely the former one. 


L = IV’XV’I, (30) 

being defined by a two-mode vector \ip) e H <g> W. For 
DP, we have the well-known Schmidt decomposition [9] 
to represent this vector, 

d d 

W) = X] = X! A "l U m v n), (31) 

i,j=l n= 1 


C. Further properties of the SEvalue equations 


Equivalent to the form of the SEvalue equations for IP 
in (20), one might formulate a second form. The solutions 
g and IF |...., bx) of the Hermitian operator L can be 
found by solving 


fl ± Lil ± \b 1 ,...,b K )=gil ± \b 1 ,...,b K ) + \x), (25) 


with the perturbation term |%), which has to fulfill for 
all j = 1,..., K and for all |a;) € H® N > an orthogonality 
relation of the form 


a;, fy+i, • • ■ ,b K \x) = 0 


(26) 


to be equivalent with the first form in (20). This treat 


ment transforms the coupled set of equations of the first 
form ( 20 ) into a single, but perturbed, eigenvalue equa¬ 


tion of the second form (251. Note that the perturbation 
|y) is an element of the (anti)symmetric subspace, since 

il ± [LU ± -gi]\b 1 ,...,b K ) = \ X ). 

Another important property of the SEvalue equations 
for IP is the behavior under certain transformations of 
the observable. Local unitaries U and shifts of L , leading 
to a transformed observable 


L' = U® N V \X 1 L + X 2 fl ± 


jj®N 


(27) 


with Ai, A 2 € R\{ 0 }, can be directly passed onto the 
solutions. If the SEvalue g together with the SEvector 
n ± |foi ,...,b K ) is a solutions of the SEvalue equations 
for IP of L, then the operator L' has the corresponding 
solutions: 


SEvalue: g' = \\g + A 2 , (28) 

SEvector: ft ± |6i,..., b' K ) = C^® jv n ± |6i, ...,b K ) 

=fl ± U® N \b 1 ,.. .,b K ). (29) 

Hence, by solving the SEvalue equations for IP for a given 
L we gain the solutions for a whole class of observables. 
Additionally, we get for Ai = 1 and A 2 = 0 that the 
SEvalues are invariant under local transformations. 


in terms of orthonormal sets {|u „)}^ =1 and {|u n )}(( =1 as 
well as non-negative coefficients X n > 0. For a Fermion or 
Boson state we get the Slater decomposition - as studied, 
for example, in Refs. naiiiiss - as 


d 

I/) =ri”IV’) = Alibi) 


i,j =1 


Ld/2j 

= ^2 K n(\w2n-l,W 2n ) - \w 2 n,W 2 n-l)), 
n—1 

d d 

|6) =f[ + |V’) = b hi\*>3) = '}2 K 'n\ w 'n, w 'n)> 

i,j =1 n—1 


(32) 

(33) 


for fij = orthonormal \w n ), denoting the 

largest integer less or equal to x, and K n > 0 ; and for 
bi,j = b jt i with orthonormal \w' n ) and n' n > 0 . 

In Appendix [C] the solution of the SEvalue equation 
for IP is explicitly computed for the considered observ¬ 
able (30). Here, let us summarize the results. It is worth 
pointing out that the nontrivial solutions, i.e., j / 0 , 
have forms which are directly related to the decomposi¬ 
tions (32) and (33). Namely, we get for Fermions 


SEvalues: g n = 2k^, (34) 

SEvectors: \w 2n -i) A | w 2n ), (35) 


and for Bosons: 

SEvalues: g n = k'2 and gk,i = ^2 + k( 2 , (36) 

SEvectors: \w' n ) V \w' n ) and \w£ t ) V \w^ t ), (37) 


respectively, with |iy^) = ■ Now, 

the entanglement condition ( 22 ), can be written in terms 
of the fidelities: 


(b\g B \b) > 


f / 2 / 2 

max < + Ki 

1 <k<l<d 1 


(38) 


or (f\g F \f) > max {^l}, (39) 

l<n<[d/2j 1 J 


for bipartite, mixed or pure entangled states of Bosons, 
Pb, or Fermions, g F . Note that, for the case of DP, we 
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also get the separable bound from the decomposition (31) 
as G = maxi< n < £ i{A^}, using the results in [52] : 


SEvalues: g n = \ 2 n , (40) 

SEvectors: \u n ) G> \v n ). (41) 


Let us apply the method to characterize the entangle¬ 
ment of the pure state which is mixed with white 
noise, 


Since the structure of (anti)symmetric product states is 
related to Bell-like states, cf. Eq. ([I]), we also consider 
Schmidt rank (SR) two states for DP. The calculation 
of the corresponding bounds is done in J5S| and applied 
in ESI- For any p in gray area of the plot SR>2, we can 
conclude that more than two tensor-product states have 
to be superimposed to describe the state (42). Thus our 
approach allows the detection of different forms of entan¬ 
glement based on a single observable. 


P = pi\ip){4>\i + (i -p) 


fi 

tr f 


(42) 


V. MULTIPARTITE EXAMPLE 


with p £ [0,1] being a noise parameter, I S {1, il + , II”}, 
and the second term being separable. For p = 0, we get 
the separable state I/tr I, which is a uniformly weighted 
mixture of all normalized product state I|ai,a 2 ) with 
1 = (ai, a. 2 |fi|ai, 02 ); cf. Eqs. ([ 2 ]) and Replacing I 

with other projections we could additionally study entan¬ 
glement for other parastatistics, e.g., for anyons El GUI, 
using the same treatment as presented for Bosons and 
Fermions. 


distinguishable, SR> 1 distinguishable, SR>2 



. 0 . 0 '. 

5 10 15 20 5 10 15 20 

d d 


indistinguishable. Bosons 



indistinguishable, Fermions 



FIG. 1: Mixing in terms of p for the state (42) depending 
on d = dim 77 is plotted. As long as p is in the gray shaded 
area, we successfully detected entanglement. The coefficients 
in Eqs. (311 and (331 for DP and Bosons, respectively, are 
chosen to be equal \k = n' k = d~ 1 I 2 (k = 1,..., d). In the case 
of Fermions, we choose Kk = (2|_d/2j) _ U 2 (k = 1,..., |_d/2j), 
see Eq. ( |32[ ), yielding a different behavior for even and odd 
dimensions d. 


In Fig. |T] we compare different quantum statistics 
regarding their entanglement properties for the mixed 
state ( |42| ) in dependence on the dimensionality of the 
single particle’s Hilbert space, d = dim"H. We apply the 
test operator in (30). As long as (L) > sup{g} (gray area 
in Fig. [T]), we have identified entanglement for the mixing 
parameter p for DP (plot: SR>1), Bosons or Fermions. 


In the following we will study a multipartite entangle¬ 
ment test, which is even independent of the spin statis¬ 
tics. We further assume dim 77 = 00 given by the or¬ 
thonormal single-mode basis {|n)}j)T 0 . The observable 
is 

L= |1,..., N)(N+ 1, ..., 2N\ + \N+1, • • •, 2N)(1, 

(43) 

For a state p, the observable L measures an interference 
term of the form 

(L)=(N+l,...,2N\p\l,...,N) 

+(l,...,N\p\N+l,...,2N). (44) 

In Appendix [D] we solve the SEvalue equations for DP 
and IP. The obtained maximal bound for A"-separable 
states is 


sup{s} = (1/2) 


K -1 


(45) 


Note that this bound is even independent of the quantum 


statistics. Therefore, the entanglement condition (12) in 


this case states: Whenever the interference (L) for N- 
particle system of DP, Bosons, or Fermions exceeds the 
bound (1/2) K ~ 1 , we have certified that the state cannot 
be K separable. 


The observable (43) may be applied to a GHZ-type 
state [60] . 


\q)=\jv{l) IX] V^WQ n \nN+l ,..., (n+l)N), (46) 


n =0 


with t £ {l,n+,n-}, ^(n^) = N !, and i/(i) = 1. 
This state is of a GHZ-type structure, because for each 
mode j holds that the individual vectors |niV + j) are 
orthonormal for different n. Using the transformations 
Tj\n)=\nN+j) in (Dl) of Appendix [73] it can be di¬ 
rectly seen that the state in (46) is a GHZ-type of state, 
| q) = fi(Ti <g) ■ ■ ■ <8 T N ) J2n=o ^ri\ n , • • • ,n). In addition, 
the pure state might be perturbed due to a randomly 
distributed parameter q ( \q\ < 1): 


P = 


'\q\<l 


d 2 qp(q)\q){q\, 


(47) 
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where p is a classical probability distribution. 

The identification of multipartite entangled Bosons 
and Fermions as well as DP is shown in Fig. [2] for a de¬ 
phasing channel, i.e., the amplitude r = |q| is fixed and 
phase ip = arg q is randomized. This uniform dephasing 
in the interval ip £ [—5, +<5] results in the density matrix 


P = 



dip 

26 


| r exp [iyj]) (r exp [i^] | 


(48) 


OO 

= ^ (1 —r 2 )r n+n sinc[<5(n—n')]^(t) 

n,n'—0 

x i|nJV+l,..., (n+l)N)(riN+l ,..., (n'+l)JV|i, 


with sinc[a;] = sin[a;]/a: (sinc[0] = 1). Hence we have for 
this state, for all I £ {1, II + , II - }, and for all K partitions 
(N u ..., Njc) the entanglement condition 

(L) = 2(1—r 2 )r sine[5] > (1/2) K ~ 1 . (49) 

As long as the expectation value in Fig. [2] is above the 
dashed lines, we certified that the state cannot be a in¬ 
separable one. Note that for a full dephasing, 5 = 7r, this 
state is diagonal in product states and, therefore, separa¬ 
ble. For no dephasing, 6 = 0, we have a pure GHZ-type 
entangled state. This example demonstrates the general 
possibility to construct spin statistics independent entan¬ 
glement tests with our approach. 



FIG. 2: Expectation value in Eq. (491 is plotted for N = 
5 (solid curve), an amplitude \q\ = l/\/3, and a uniformly 
distributed phase in the interval arg q £ [— 5, 5]. If (L) is above 
the dashed line K, then the state cannot be a A'-separable one 
- independent of the quantum statistics. 


VI. CONCLUSIONS 

In summary, we derived a method which allows the 
construction of entanglement probes in systems of Bosons 
and Fermions. These necessary and sufficient conditions 
are capable of determining full and partial entanglement 


for any number of particles. The optimization of these 
criteria is based on a set of generalized eigenvalue equa¬ 
tions. These equations yield a structural unification of 
entanglement for distinguishable particles, Fermions and 
Bosons, and it can be generalized to other parastatistics. 

We analyzed and compared a number of examples for 
studying the differences and similarities between full and 
partial separability as well as entanglement in systems of 
distinguishable and indistinguishable particles. For in¬ 
stance, the determination of entanglement in bipartite 
and multipartite as well as discrete and continuous vari¬ 
able quantum systems demonstrate the wide range of ap¬ 
plications of our technique even in the presence of noise. 
This also shows that our method is not limited to bi¬ 
partitions or small numbers of particles. Moreover, the 
construction of spin-statistics independent entanglement 
probes has been established. Such witnesses can detect 
entanglement of a quantum system independent of sym- 
metrization effects. 

We presented an approach which allows one to con¬ 
struct entanglement criteria, in principle, from almost 
all observables. However, the question which measurable 
quantity is able to witness the entanglement of a partic¬ 
ular state is open and requires further studies. Beyond 
the here presented full analytical approach, numerical im¬ 
plementations may allow one to generate more sophisti¬ 
cated entanglement probes for detecting entanglement in 
more general scenarios. Therefore, we believe that our 
approach will provide a versatile tool to characterize en¬ 
tanglement in future experiments, with applications to 
Bose-Einstein condensates or ultra-cold Fermi systems. 
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Appendix A: Symmetrization of operators 


For an IV-fold Hilbert space 'H ® N , the projection op¬ 
erators II + and 11“ are defined as 


n* 


E 

<j€lSn 



(Al) 


with P a \ai,a N ) = \a a(tl ),..., a CT (jv)) for any permu¬ 
tation a £ Sn■ It holds (IF^t = fl*. We may study 
Hermitian operators in a product basis operator expan¬ 
sion, given by terms of the form 

X = Yi (g> • • • ® Yn, (A2) 

with Yj = Yj for j = 1,..., TV. The symmetric form of 
X is defined as 

1(sym) = E ® •• • ® t(iv). (A3) 
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We claim 

fr^Yn* = i-( s y m )n± = n±jY (sym ). (A4) 

The first equality can be directly computed, since for all 
|ai,..., a_/v) holds: 


N 


N 


n ± xn ± 

i =i 

(±l)l cr l + l T I 


> = n± v MA 

' ^ N\ 

o-eSjv j =i 

N N 


hl a o-(j)) 


= E 

c t,t£Sn 


= y — 

^ N\ 

res N j =i 


N\ 2 

N 


3 = 1 


Y, 


O') 


r(<r(j))) 


3 = 1 


W £ 


AT 


H^zSn 


N\ 


3 =1 


AT 


=J £(sym)n± (O) | a v 


(A5) 


f=i 


where we used a substitution /r = r o a and (±i)l r l+l cr l = 
(±l)l roo 'L The second equality in (A4| follows from the 


fact that X and EF are Hermitian operators, 

n±^( S ym) = n±xn± = (ri ± xn ± ) t 
=(n ± i' (sym) ) t = 


(A6) 


For X = i 


I, we get from (A4) that II ± is 


idempotent. An observable L which solely acts on the 
corresponding subspaces ("H viv or T-L aN ) should fulfill the 
commutation relation [L, IF] = 0. From (A4) follows 
that this is fulfilled for every L = l/ sym ). 


Appendix B: Existence of partial separable Bosons 
and Fermions 


with the orthogonal vectors |$o) = |1, 2) ± 12,1) + |3,4) ± 
|4,3), |$!) = 12, 0) ± 10, 2), |$ 2 > = |0,1) ± |1, 0), |<i> 3 ) = 
|4,0) ± |0,4), and |<J> 4 ) = |0,3) ± |3, 0). The reduced state 
- tracing with respect to the first subsystem - is 

/5 re d = tri|4' ± )(T ± | (B4) 

J$o)($o| + |$l)($l| + |$ 2 )($ 2 | + |$3)($3| + |$4)($4| 

6 

and it has a rank of five, rank(p re d) = 5. 

For proving that the states | ) cannot be fully sepa¬ 

rable, IIP) £ K)V|a 2 )V|a 3 ) and |T~) ^ |ai)A|a 2 )A|a 3 ), 
let us study the properties of the reduced density matrix 
of fully separable Boson and Fermion states. We have 

|s ± ) =n ± (|ai) <g> |a 2 ) <g) |a 3 )) 

= ^K) <S> |si) + ^|a 2 ) <8 |s 2 ) + ^|a3) ® |-S3>, (B5) 

with |si) = |a 2 ,a 3 ) ± |a 3 ,a 2 ), |s 2 ) = |a 3 ,ai) ± |ai,a 3 ), 
and |s 3 ) = |ai,a 2 ) ± |a 2 ,ai). Consequently, the range or 
image of the partially reduced operator is 

7?.=Im(tri|s ± )(s ± |) C span{|si), |s 2 ), \s 3 )}=H' , (B6) 

where the linear span TV has a dimensionality of three 
(if |si), |s 2 ), |s 3 ) are linearly independent) or less. Since 
dim 1Z < 3, it follows that the rank of the partially 
reduced operator of any fully separable state of IP is 
bounded by three: 

rank(p(. ed ) < 3, for all (B7) 

Pred= tr i n ± |ai,a 2 ,a 3 )(ai,a 2 ,a 3 |n ± . 

Finally, we conclude that the states |4 f± ) cannot be fully 
se parab le because of rank(p re d) = 5 > 3. Thus, the states 
in (B1), using Eq. are authentic examples of partially 
separable states of Fermions and Bosons that are not fully 
separable. 


Let us prove that the set of partially separable states 
of IP includes more elements than the fully separable 
ones. For this reason, we consider the orthonormal basis 
{|0),..., |4)} and the three-partite vector In') G (C 5 ) ® 3 
in Eq. (J9|) . It is partially separable in the tensor product, 
i.e., |’P) = |0) <B> |$). Applying the symmetrization or 
anti-symmetrization operator, II , we get 

|T+)=n+|4')G(C 5 ) V3 and |4'-)=II"|4 ')g(C 5 ) A3 , (Bl) 
or, more explicitly, we have the expansion 

|T±) =|[|0,1,2) ± |0,2,1) + |0, 3,4) ± |0,4, 3) 

+ |1,2,0)±|1,0,2) + |3,4,0) ±|3,0,4) 

+ |2,0,1) ± 12,1,0) + |4,0,3) ± |4,3,0)] (B2) 


-x[|0) ® |$ 0 ) + |1) ® |$i) + |2) 8 |$ 2 ) 

0 

+ |3) ® |$ 3 ) + 14) 8 |$ 4 >], 


Appendix C: Bipartite observable 

In the following two subsections, the solution of the 
SEvalue equation for IP will be explicitly computed for 
the considered observable (30). For simplicity, we will as¬ 


(B3) 


sume in the following L = |/)(/| (L = |6)(6|) for the two- 
Fermion (two-Boson) system which yields L = II - LII - 

(Z = n+Ln+). 


1. Solution for Fermions 

First, we give a unitary state representation for arbi¬ 
trary d-dimensional (d € N U {oo}) pure states of two 
Fermions. We start with a Fermion state, 

d 

I/) = E /bilbi). with fi,3 = 

i,3 =1 


(Cl) 

















and introduce the skew-symmetric coefficient matrix 

M f = (fijhj = ~Mj. (C2) 

Using the Autonne-Takagi factorization in Ref. [55], we 
find the Slater decomposition of this coefficient matrix, 


M f = UDU T , with ffu = i 


(C3) 


d/2 


and (d even): D = © Kj 
i=i 

(d- 1)/2 

or (d odd): D = 0 * 

3 = 1 


0 +1 
-1 0 

' 0 +1 
-1 0 


®[°]i 


with D being a block diagonal matrix containing anti¬ 
diagonal 2x2 blocks and Kj > 0. This yields Eq. (32) in 
the form 


Ld/2j 

|/) = U <g) U Y «n(|2 n - 1, 2 n) - |2n, 2 n - 1)), (C4) 

n—1 

with an orthonormal single-mode basis {11), |2),...} and 
U\k) = | Wk)- Since the SEvalues are invariant under 
unitary separable operations U , we assume, without 
loss of generality, that U = 1. 

Now, we consider more general operators having an 
expansion as 

Ld/2j 

L = Y^ A TOi „(|2to — 1,2m) — |2m, 2m — 1)) 

m,n =1 

x ((2n-l,2n| - (2n,2n-1|), (C5) 


which includes the special case L m ,n = K m K n of pro¬ 
jection operator L = |/)(/|. Since L = If - All - , the 
SEvalue equations for Fermions in the second form (25) 
read as 

L|ai,a 2 ) = S^(|ai,a 2 ) - |a 2 ,ai)) + |y). (C 6 ) 


Using 7 „ = ((2 n — 1, 2n\ — (2 n, 2 n — l|)|ai, a 2 ), we get 


for all n = l,...,[d/2j. For the special case L m ^ n = 
KmKm we get the maximal SEvalue as 

G = max{2 kV\. (C9) 

n 

Note that trivial solutions, <7 = 0, can be obtained by 
\k) A 1 1) with (k, l) 7 ^ (2 n, 2 n + 1). 


2. Solution for Bosons 


Again, we first give the state representation according 
to Ref. )55j for arbitrary d-dimensional pure state of two 
Bosons. This means that the symmetric state 

d 

I b) = Y Kj with bij = bjj, (CIO) 

i,j =1 


can be identified with a symmetric coefficient matrix 

Mb = (bij)ij = Mb and M b = UDU T , (Cll) 


with U^U = i and D = diagfc :),... ,n' d ] > 0. Thus the 
symmetric Slater representation (33) of the state is 


d 

| b) = U K ' n \ n,n), with U\n) = \ w' n ). (C 12 ) 

n=1 


As in the previous example for Fermions, let us con¬ 
sider the more general operator 


d 

L= Y L m ,n\m,m)(n,n\. (C13) 


m,n—1 

The relation L = n + LII + simplifies the SEvalue equa¬ 
tions for Bosons in the second form to 

L\a 1 ,a 2 ) = fl^(|ai,a 2 ) + |a 2 ,ai)) + |x)- (C14) 

Using 7 „ = (n, n|ai, a 2 ), we can now write 
d f d 


L\a 1 ,a 2 ) = Y 


m= 1 




n—1 


(C15) 


L\a\, o 2 ) 


Ld/2J 

= £ 


m—1 


Ld/2J 

^ ) A m 7l 7n 

n=l 


(|2 m 


(C7) 


1 , 2 m) — | 2 m, 2 m — 1 )). 


We fi nd t hat L|ai,a 2 ) is already diagonalized in the 
form ( |C4[ ). Hence, the orthogonality of |ai,a 2 ) to the 
perturbation |x) is fulfilled if 

ft” |a-i, a 2 ) =—(|2n—1,2n) — \2n, 2n— l))=|2n— l)A|2n), 
g =2 Tn,ra, (C8) 

lx) = Y L m,n{\ 2m - !> 2m) - 12m, 2m - 1)), 

m^n 


Hence one class of solutions with |ai) = |a 2 ) is given by 

n + |ai,a 2 ) =| n,n) = |n) V |n), 

<? =£n,n, (C16) 

lx) = Y L m,n\m,m). 

m^n 


Unlike in the Fermion case, we have to take a brief 
look on the decomposition of product states of Bosons. 
Namely the state If + |a i,a 2 ) for any |ai) ^ |a 2 ) has a 
decomposition, cf. Eq. (|C12), as 


n+| ai , o 2 ) =U’ <g> U'( A;|l, 1) + A' |2,2», 
for |o 1(2) ) =& [ v / A[|l)+(-)ivE|2)' • 


(C17) 
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Hence, we get a more involved set of solutions of 
Eq. (C15) in the form (for k ^ l): 

ri+|ai,a 2 ) =X' k \k, k) + X[\l, l), (C18) 

lx) = ( L m,kX' k + L m,iX’ l )\m,m), (C19) 

m^k,l 

where the coefficients A' fe a nd X[ have to be determined. 
We insert (C18) and ( |C19 ) into ( |C14[ ), 


Lri + |ai,a 2 ) - |x) = gn + |ai,a 2 ), 


(C20) 


with z^(II ± ) = N\ and z/(l) = 1. Due to this fact, we may 
define the AT-separable vectors 


|ui,...,UK)=|l,...,A r ) 

and \w-l, ... ,w K )=\N+l, .. . ,2 N), 


(D4) 


which are orthogonal for Fermions, Bosons, and DP and 
any partition (Aq,..., Nk)- 

As a last fact before we solve the SEvalue equations 
for this operator, let us recall an example of the standard 
eigenvalue problem: 


and find that the remaining terms to be computed are 
(Lk,kX k + L k iX[)\k, k) 

+(Li,kX' k + L lt iX't)\l, l) = g(X k \k, k) + X[\l, l)). (C 21 ) 

This is a standard eigenvalue problem in C 2 , which has 
the solutions 

9 ± = Lfc ’ fc +^ M±A , A = ^(L fe , fc -L M ) 2 +4|L fcji p, 

X' k =2L kyl and A) = L M - L Kk ± A, (C22) 

with the Hermiticity condition Iq tk = L k j. 

Again, in the particular case L m ^ n = we get the 

simplified solutions 

9 = g~ = 0, and g + = n' k + nf. (C23) 

Combining the solutions of the form |ai) V |ai) and |oi) V 
|a 2 ) as well as using the fact that n' k + k[ > k[ , we get 
the maximal SEvalue as 

G = max{A^ + A 2 }. (C24) 


M =m\m w ){m v \ + m*\m v ) (m w \, (D5) 

|M±) = (\m w ) ± ^j- |m„)^ /V2, (D 6 ) 

H±=±\m\, (D7) 

with complex to 7 ^ 0 , orthonormal {\m w )^ |to„)}, and 
|/x±) being the eigenvectors of M to the eigenvalues p±. 

Now, let us use the first form of the SEvalue equation 
for IP and DP of the operator ( |43| . Since the spanned 
subspace of L is span{|vi,..., vk), |uq, • ■ •, wr-)}, let us 
expand 


\bj) = Pv,j\vj) + 0w,j\wj)- (D 8 ) 


We get for the jth SEvalue equation the two components 


= 5 (^ 1 ( 1 ) 571 ^), 

(wjKILI)^) =g(w 3 \(l)^\b J ). 

Equivalently, we get by a rescaling with i/(I) 


(D9) 


Appendix D: Multipartite observable 


We considered an interference operator L in Eq. (43) 


whose expectation value is the real part of an off- 
diagonal element of the density operator p, ( L ) = 
2 Re(/ 9 ( 1> ... ) jv),(jv+i...., 2 JV))- Local unitary operations al¬ 
low the generalization to other off-diagonal elements or 
phase shifts, Re{exp[iip]P(i,...,N),(N+i,...,2N))- 

The injective transformations T) of the orthonormal 
basis {|n)}„ eN , 


fj\n)=\nN+j) for j = 1 ,..., N, 


(Dl) 


n (pv,iPw,i) p w ,j=g n d^yi 2 +i^.ii 2 ) Pv,j, 

i¥=J *¥=3 

n (Pw,iPv,i) pv,j=g n (\Pv,i\ 2 +\PwA 2 ) Pw,j, 


(DIO) 


which has the structure of the eigenvalue problem in (D5) 


with the solution in Eqs. (D 61 and (D7). Hence for each j 


we get the solution for components with \/3 v j\ = \Pw,j\ = 
1 /x/ 2 , yielding |/3„,i | 2 + \Pw,i\ 2 = 1 anc l the eigenvalues 


g = ± \H(p* W ' i p v , i )\=±(l/2) 


K -1 


(Dll) 


are constructed such that one can directly see that the 
for all j.j' = 1,..., N and n, n' G N an orthogonality is 
given, (nN + j\n'N + j') = Therefore, we get 

for I G {l,n+,II-} the orthogonality relation 

(1,...,N\I\N+1,...,2N) =0 (D2) 

as well as the normalizations 

(N+l,...,2N\I\N+l,...,2N) 
={l,...,N\I\l,...,N)=l/v$), (D3) 


Note for ]([,.,■ = 0, we get the trivial SEvalue 

g = 0 and, for example, the SEvector S|6 i,..., bx) = 
%i,..., vk)- Finally, the maximal SEvalue is 

G = sup{ 5 } = (1 /2 ) if " 1 . (D12) 

Note that this result is independent of the particular Re¬ 
partition (A r i,... , Nk ) and, due to especially chosen or¬ 
thonormality in ( |D2[ ), the result is also independent of 
the spin statistics. 
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